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ABSTRACT. Linear programming bounds provide an elegant method to prove 
optimality and uniqueness of an (n, N, t) spherical code. However, this method 
does not apply to the parameters (4, 10, 1 /6). We use semidefinite programming 
bounds instead to show that the Petersen code, which consists of the midpoints 
of the edges of the regular simplex in dimension 4, is the unique (4, 10, 1/6) 
spherical code. 



1. Introduction 

Let C be an iV-element subset of the unit sphere S 1 ™ -1 C R n . It is called an 
(n, N, t) spherical code if every two distinct points (c, d) of C have inner product 
c • d at most t. An (n, N, t) spherical code is called optimal if there is no (n, N, t') 
spherical code with t' < t. 

Only for a few parameters optimal spherical codes are known. The table O 
page 115] lists all known cases in dimension n = 3. The tables [ 17 , Table 9. 1] and 
iPTTl Table 1] list all known cases in which optimality can be proven using linear 
programming bounds. 

One source of optimal spherical codes are iterated kissing configurations com- 
ing from the Eg root lattice in dimension 8 and the Leech lattice in dimension 24 
(see H3). Starting from the sphere packing defined by these lattices one fixes one 
sphere and considers all spheres in the packing touching the fixed one. The touch- 
ing points, also called a kissing configuration, form (8, 240, 1/2) and respectively 
(24, 196560, 1/2) spherical codes. Then one views the kissing configuration as a 
packing in spherical geometry and repeats this construction. One gets (7, 56, 1/3) 
and respectively (23, 4600, 1/3) spherical codes. 

More formally, one picks a point x S C from an (n, M, 1/k) spherical code 
C in which x has M' points N x C C with inner product \/k. Then the points 
(N x - x/k) /a/1 - 1/k 2 form an (n - 1, M' , l/(k + 1)) spherical code. 
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In this way one gets sequences of spherical codes with parameters 

(8, 240, 1/2), (7, 56, 1/3), (6, 27, 1/4), (5, 16, 1/5), (4, 10, 1/6), (3, 6, 1/7), 
and 

(24, 196560, 1/2), (23, 4600, 1/3), (22, 891, 1/4), (21, 336, 1/5), (20, 170, 1/6). 

By using linear programming bounds Levenshtein ifTTl proved that every sharp 
(see Section [3]> spherical code is optimal. Levenshtein's theorem applies to all 
spherical codes above except to those with parameters (4, 10, 1/6), (3, 10, 1/7), 
(21, 336, 1/5), (20, 170, 1/6). In all optimal cases the spherical code is also unique 
up to orthogonal transformations. This was proved for the cases (8,240,1/2), 
(7, 56, 1/3), (24, 196560, 1/2), (23, 4600, 1/3) by Bannai and Sloane O and for 
(22, 891, 1/4) by Cuypers HH and independently by Cohn and Kumar H3 (who 
also corrected a minor error in the (23,4600,1/3) case). For the cases (6,27,1/4), 
(5, 16, 1/5) see the discussion in IfTTl Appendix A]. One should point out that 
optimality does not imply uniqueness as one can see from the sharp (q(q 3 + l)/(q+ 
1), (q + l)(q 3 + 1), l/g 2 ) spherical codes from [10]. For some q there are two 
different spherical codes with these parameters. 

Based on massive computer experiments Cohn et al. JU Section 3.4] conjectured 
that the (4, 10, 1/6) spherical code is optimal and unique. As we explain in Sec- 
tion [2] this spherical code is closely related to the Petersen graph and we call it the 
Petersen code. Whether the above spherical codes with parameters (21, 336, 1/5) 
and (20, 170, 1/6) are optimal and unique is currently unclear. At least in all these 
cases linear programming bounds cannot be used to show optimality. A (3, 6, 1/7) 
spherical code is not optimal because the vertices of the regular octahedron form a 
(3, 6, 0) spherical code which is a sharp spherical code. 

The main result of this paper is the following theorem which proves the conjec- 
ture. 

Theorem 1.1. The Petersen code is an optimal (4, 10, 1/6) spherical code. Up to 
orthogonal transformations it is the unique spherical code with these parameters. 

The proof is based on the semidefinite programming bounds for spherical codes 
developed in and 0. Currently this is the only new case we know where the 
semidefinite programming bound is tight and the linear programming bound is not. 
Another known case seems to be 8 points in S 2 which was solved by Schiitte and 
van der Waerden in fl9l . The linear programming bound gives 8.29 whereas our 
numerical calculations suggest that the semidefinite programming bound is tight. 

We could not prove optimality of (21,336,1/5) and (20,170,1/6) spherical 
codes using semidefinite programming bounds. For the first case the linear pro- 
gramming bound equals 392 whereas our numerical calculations suggest that the 
semidefinite programming bound is approximately 363. However, we run into se- 
rious numerical problems here and at the moment we cannot definitely rule out 
that the semidefinite programming bound is sharp. For the second case the linear 
programmig bound and the semidefinite programming bound coincide: They both 
give 206.25. 
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The structure of the paper is as follows: After giving some constructions and 
properties of the Petersen code in Section |2l which also reveal the origin of its 
name, we show in Section [3] that one cannot prove Theorem 11.11 using linear pro- 
gramming bounds. In Section 0] we recall the semidefinite programming bounds 
and in Section[5]we present a proof of Theorem 11.1 1 based on them. 

2. Constructions and properties of the Petersen code 

There are many possibilities to construct the Petersen code and we already gave 
one. Here we give two more. 

The next construction justifies the name "Petersen code". The Petersen graph is 
a graph with 10 vertices and 15 edges. The vertices are given by the 2-element sub- 
sets of a 5-element set and they are adjacent whenever the corresponding 2-element 
subsets have empty intersection. Every point of the Petersen code corresponds to 
a vertex of the Petersen graph and the inner product between two points is —2/3 
whenever the corresponding vertices are adjacent. The inner product is 1/6 when- 
ever the corresponding vertices are not adjacent. This defines a Gram matrix having 
rank 4 which is unique up to simultaneous permutation of rows and columns. The 
number of ordered pairs in the Petersen code with inner product —2/3 is 30 and 
those with inner product 1/6 equals 60. 

In the Petersen graph every vertex has three neighbors, every pair of adjacent 
vertices has no common neighbors and every pair of nonadjacent vertices has ex- 
actly one common neighbor. So it is a strongly regular graph with parameters 
v = 10, k = 3, A = 0, \i = 1. It is easy to see that it is uniquely defined by these 
parameters. For more information about strongly regular graphs see [4] and [8 ]. 

The next construction is geometric: After applying a suitable similarity transfor- 
mation the midpoints of the edges of the regular simplex in dimension 4 form the 
Petersen code. Sometimes, this construction is also called the second hypersimplex 
A(2,5). The second hypersimplex is the 4-dimensional poly tope defined as the 
convex hull of the points e« + e 7 with 1 < i < j < 5 where is the z-th standard 
unit vector in M 5 . For more information about second hypersimplices see lfl~8l . 

By lfl~5l Theorem 5.5] the Petersen code forms a spherical 2-design: A spherical 
code C C S n ~ 1 forms a spherical M-design if for every polynomial function 
/ : W 1 — > R of degree at most M, the average over C equals the average over the 
sphere S n ~ l . 

3. Linear programming bounds 

Linear programming bounds provide an elegant method to prove optimality and 
uniqueness of an (n, N, t) spherical code. In particular a theorem of Levenshtein 
lfi"7l Theorem 1.2], which covers many cases in a unified way is based on them. 
Before we prove that linear programming bounds cannot prove the optimality of 
the (4, 10, 1/6) spherical code we briefly review the underlying notions (see also 
e.g. [15, Theorem 4.3], 021, 03 Chapter 9], E Theorem 2.1]). 

The positivity property of the Gegenbauer polynomials C^ 2 1 (see 01 Chapter 
6.4]), which are normalized by C^ 2 = 1, underlies the linear programming 
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bounds for spherical codes in S n : For every degree k = 0, 1, . . . and every finite 
subset C of 5 n_1 we have 

(1) £ c^V <0>o. 

(c,c')ec 2 

One formulation of the linear programming bounds is as follows. 

Theorem 3.1. Let F(x) be a polynomial with expansion 

d 

(2) F{x)=Y J fkC n k l2 -\x) 

k=0 

in terms of Gegenbauer polynomials C^ 2 \ Suppose that 

(a) all coefficients fj~ arc nonnegative, 

(b) /o > 0, 

(c) F(x) < for all x G [-l,t]. 
Then an (n, AT", t) spherical code satisfies 

(3) N < Fm 



A 

Proof. For an (n, A 7 ', i) spherical code C we have the inequalities 

(4) NF(1)> E F ( c " c ') = E ^(c-c')>iV 2 /o. 

(c,c)ec 2 (c, c ')eo2 (c,c')eC 2 

where the first inequality is due to (c) and the second due to (a) and the positivity 
property (GQ). This together with (b) implies ®. □ 

If there exists an (n, N, t) spherical code C so that N = L^(l)//oJ in ©> then, 
of course, C is a maximal (n, N, t) spherical code, i.e. N is the maximal number of 
points one can place on the sphere S™' 1 so that distinct points have inner product 
at most t. If furthermore = F(l)//o, then C is an optimal (n, N, t) spherical 
code. This can be seen as follows. If ([3]) is tight it follows from the proof that for 
an (n, N, t) spherical code C one has F(c ■ d) = for distinct c, d G C. Suppose 
C' is an (n, N, t') spherical code with t' < t. Then, F(c ■ d) = for all distinct 
c, c' G C . Now we perturb C continuously to another (n, N, t") spherical code 
C" with t' < t" < t. Still we would have that c • d is a root of the polynomial F 
for all distinct c, d G C" yielding a contradiction. 

Levenshtein's theorem says that for every sharp spherical code there is a polyno- 
mial satisfying the assumptions of Theorem 13.11 for which ([3]> is tight. A spherical 
code C is called sharp if it is a spherical M-design and the number m of different 
inner products between distinct points satisfies M > 2m — 1. The Petersen code 
is a spherical 2-design and there are 2 different inner products between distinct 
points. Thus, Levenshtein's theorem does not apply to it. 

Now we show that it is not possible to prove the optimality of the Petersen code 
with help of Theorem l3 . 1 1 Suppose that the polynomial F(x) = l+Ylt=i fkC\{x) 
satisfies f k > for k = 1, . . . , d and F(x) < for all x G [-1, 1/6]. If F would 
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prove that the Petersen code is optimal, then the inequalities in © are equalities, 
so we would have that 

d 

(5) 10 = F(1) = 1 + 

k=l 

and that 

(6) = F(-2/3) = F(l/6), 
and furthermore that for all fc with /fc > 

(7) 0= £ C fc 1 (c-c , ) = 10 + 30C fc 1 (-2/3) +60^(1/6). 

(c,e')eC 2 

We shall show that © only holds for k = 1 and = 2: By fl] (6.4.1 1)] we have 
the following expression 

k 

(8) C£(co80) = _^£cos((fc-2i)0). 

i=o 

Hence, 

(9) lim C*(-2/3) = lim (7^(1/6) = 0, 

k— >oo fe^oo 

so that for sufficiently large /c, ([7]) cannot hold true. Checking the remaining cases 
it follows that Q is only valid for A: = 1,2. Hence, F is of degree 2, but then F 
cannot satisfy the conditions (f5]) and (© and F(x) < for x € [—1, 1/6]. 

This argument gives rather pessimistic estimates. In fact, numerical computa- 
tions suggest that for all d > 3 the optimal polynomial is 

„, . 2270 2775 ,4 2 1, 1500 , o N 

(10) = H xH (—x ) H (2a; 3 -x), 

v ' v ; 680 680 v 3 3 ; 680 v ; 

and so the best upper bound one can probably prove using Theorem |3.1| is 10.625. 

We checked this for all d < 40 by computer. 

4. Semidefinite programming bounds 

As we have seen above the positivity property of the polynomials C^ 2 1 plays a 
crucial role for the linear programming bounds. For the semidefinite programming 
bounds this is replaced by the positivity property of the matrices S%. From we 
recall the matrices Su and their positivity property. First we define the entry (i, j) 
with i, j > of the (infinite) matrix YJ} containing polynomials in x, y, z by 

{Y k n ) hj (x,y,z)=xY ■ 
(H) ((l-^Xl-y 2 ))^/ 2 - 3 / 2 ( ^ ^ 



V(l-x 2 )(l-y 2 ) 



and then we get S% by symmetrization: 

(12) S3f = ^E ff ^ 
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where a runs through all permutations of the variables x,y,z which acts on the 
matrix coefficients in the obvious way. The matrices S% satisfy the positivity prop- 
erty: 

(13) for all finite C C S n '\ ^ S%(c ■ J, c ■ c", c' ■ c") y 0, 

(c,e',c")eC 3 

where u y 0" stands for "is positive semidefinite" where we mean that every fi- 
nite minor is positive semidefinite. Note that the difference between (fTTT> and the 
original El (11)] is due to a change of basis which does not affect the positivity 
property. 

The interval [— l,t] of the linear programming bounds is supplemented by the 
domain 

(14) D = {(x, y, z) : -1 < x,y, z < t, 1 + 2xyz - x 2 - y 2 - z 2 > 0}. 

We need some more notation. The space of (finite) symmetric matrices is a 
Euclidean space with inner product (F,G) = trace (FG). The cone of positive 
semidefinite matrices is self dual, i.e. one has (F, G) > for all positive semidefi- 
nite G if and only if F is positive semidefinite. If F is a symmetric matrix with m 
rows and m columns, then we interpret (F, as the inner product of F with the 
principal minor of S% of appropriate size. 

Now we can state the semidefinite programming bounds. The following poly- 
nomial formulation can be deduced from \2, Theorem 4.2]. We provide an inde- 
pendent proof which has the additional feature that it gives information in the case 
when the theorem provides tight results. 

Theorem 4.1. Let F(x, y, z) be a symmetric polynomial with expansion 

d 

(15) F(x,y,z) = J2( F k,S%), 

k=0 

in terms of the matrices S"^. Suppose that 

(a) all Fk are positive semidefinite 

(b) Fq — /qEq y for some fo > (Eq is the matrix whose only nonzero 
entry is the top left corner which contains 1 ), 

(c) F(x, y, z) < Ofor all (x, y, z) £ D, 

(d) F(x,x,l) <Bforallx£ [-l,t]. 
Then an (n, N, t) spherical code satisfies 



(16) N < — . 

Proof. Let C be an (n, N, t) spherical code. Define 

(17) S= F(c-c',c-c",c' -c"). 

(c,c',c")GC 3 

Split this sum into three parts according to the indices C\ , C2 , C3 C C 3 where 
Ci contains all triples with i pairwise different elements. The contribution of C\ to 
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S is NF(1, 1, 1), the one of C 2 at most 3N(N - 1)B and the one of C 3 is at most 
zero. Together, 

(18) S < NF(1, 1, 1) + 3N(N - 1)B. 
On the other hand, 

d 

(19) S = J2(F k , SUc-c',c-c",c'-c") 

(20) > {f Eo, Yl SUc-c',c-c",c' -c")) 

(c,c',c")GC* 3 

(21) = iV 3 /o, 

yielding the statement of the theorem. □ 
A few remarks about the theorem and its proof are in order. 

If the bound (fl6l ) is tight, then all inequalities in the proof must be equalities. In 
particular, the univariate polynomial F(x, x, 1) — B has roots at the inner products 
c • d for distinct c, d G C. So we can argue in the same way as in the case of the 
linear programming bounds that tightness implies optimality. 

If the bound ( fl6l ) is tight, we have the following identities: Let C be an (n, N, t) 
spherical code with 

D{C) = {(c • d, c ■ c", d ■ d') : (c, d, c") G C7 3 }, 

U-) 1(C) = {c-d : (c,d) G C 2 ,c^d}. 

Let F be a polynomial satisfying the hypothesis of Theorem 14. II with constants B 
and f and proving the tight bound (SB + ^9B 2 + 4/ (F(l, 1, 1) - 3B))/2/ . 
Then 

(i) Ar 2 /o-^(l,l,l)-3(^-l)B = 0, 

(ii) y, z) = for all (x, y, z) G 73(C), 
(hi) F(x, x, 1) = B for all a; G 1(C), 

(iv) (i^A, E( C C ' c'')ec3 ^(c • c', c • c", c' • c")) = for all A: = 1, . . . , d, 

(v) (F , E(c,J, c »)eC3 Sq(c • c', c • d', d ■ d')) = N 3 f . 

Semidefinite programming bounds are at least as strong as linear programming 

bounds: If G = Ylt=o 9k^k ^ s a P°ly nom i a l which satisfies the hypothesis 

of Theorem l3.1[ then the polynomial F(x, y, z) = (G(x) + G(y) + G(z))/3 satis- 
fies the hypothesis of Theorem 14. II with B = G(l)/3 and /o = 5o- This is because 
one sets Fo = goEo and from (2J Proposition 3.5] it follows that one can express 
G with semidefinite matrix coefficients. 

From [3 . Lemma 4.1] it follows that one can express every symmetric polyno- 
mial in the form (fT5l) . However, this expansion is not unique, e.g. 

(23) x + y + z=(( 3 ° /2 f),SZ) + ((o),S?) 

= (( ° 3 ),sz) + {(3),s?), 
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where only the second expansion involves semidefinite matrices and where 

( 24 ) ss = ( {x+ ; +z)/3 ( j;;r + i / ) 3 /3 ) > s i = ( (^.+^)/3-(w + ^)/3 ) 



5. Proof of optimality and uniqueness 

In this section we prove Theorem 11.11 with the help of Theorem 14.11 Although 
we can present a proof which one can verify essentially without using computer 
we relied heavily on computer assistance to find it. 

To show that the Petersen code is the unique (4, 10, 1/6) spherical code we use 



the matrices Fq € 



p4x4 



(25) 



F n 



Fx 



Fi € M ax3 , F 2 € 

/2882/3 114 - 

114 324 

2500 216 





3588 
-4536 



Dlxl 

-2500 
216 
8716 
1296 



given by 

\ 


1296 
11664/ 




,F 2 



(2000) . 



Let 

(26) 



m 



i5X 



< i < j < k, 



where a runs through all permutations of the variables x, y, z, be the polynomial 
which one gets by symmetrizing x l y J z k . Then, 



F(x, y, z) =11664m 32 o + 11664m 22 i + 7128m 220 - 9072m 2 n 

+ 432m 2 i - 2412m m + 324mu + 228mi 00 - 118/3, 



(27) 
and 

(28) F(x,x,l 



B 



1 



3888 



2 \ 2 / 
x + -J ( :r 



1 



20 



l)\ X + 9 X+ 27 



It is a straight forward computation that F satisfies the condition of Theorem 14.11 
with F{1, 1, 1) = 59750/3, B = 250, f = 800/3 so that it shows N < 10 for a 
(4, N, 1/6) spherical code. This finishes the proof of the optimality. 

Before showing uniqueness, let us describe how we derived Fq, F\, F 2 . We have 



(29) 



So(x,y,z) 

Sf(x,y,z) 
S$(x,y,z) 



( 1 

"T-100 

"T-200 "1210 

\m 30 o "1310 

("T-100 - "1110 
m 110 - m 2 i 
m 2 io 



TtlOO "1200 "i30o\ 

"1110 "1210 "1310 

m 2 20 "1320 

m 32 "1330 / 

"1-110 

m m 

"1211 - "1320 



"1210 "1210 - "l31o\ 
"1220 "1211 - "1320 , 



"1310 



"1221 - "1330 / 



-\ + |m 20 o - 3mm + "I220 ■) 
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ThenO = ELo(^.^) for 



(30) K 



1,0 



f° 


1 

2 







1 

2 


1 




























0/ 



1,1 




(0) 



K 



2,0 



K 



3,0 



K, 



4,0 



( \ 







-! °\ 







I 2 
\ 



-1 







.#2,1 




2.2 



/o 




Vo 
/o 






0\ 

1 




K. 



3,1 



1 



, K; 



3,2 



(o) 



0/ 






0; 



0\ 

-i I 

2 2 






4,1 



0/ 



1 r 

I . K 

1 0, 



4,2 



(o) 



i.e. the matrices form a basis of the kernel of the linear map which assigns 
symmetric polynomials to the matrix coefficients. From the discussion following 
the proof of Theorem 14 . 1 1 we know that the matrix entries have to satisfy the equal- 
ities (i)-(v) where 



/1000 250 



S*(c-c',c-c",c' -c") 



(c,c',c")&C 3 



(31) 



b 1 {C ■ C ,C ■ C ,C ■ C J 



(c,c',c")&C 3 




125 \ 
9 \ 


125 

M, 

648/ 



^2 Sf (c • c', c • c", c ■ c") = (0) . 

(c,c',c")&C 3 

We restrict our search to polynomials F satisfying 
dF / 2 2 h OF/ 2 11 



(32) dx( :Y :Y() 



( 

' dx V 3 ' 6 ' 6 



dF / 2 2 



0. 



Furthermore, we restrict our search to those polynomials lying in the subspace of 
dimension 9 spanned by 



(33) 



W320 , ^221 , ™220 , ™211 , ™210 , "^111 , »™110 , ™100 , 1 ■ 
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The one dimensional affine subspace 

Ey(x, y, z) = (ll664m 32 o + 9720m 220 - 1296m 2 i - 6480mm 

+ 2268mno - 108toi 00 - 18) + 7(34992m 22 i - 7776m 220 

( ' 4 ' - 27216m 2 n + 5184m 2 io + 12204mm - 5832mno 

+ 1008mioo - 64) , 7 el, 

satisfies all these linear equalities. We have 

2 

F-y(x,y,z) = ^2(A k , 



(35) 



,St) +1 (B k ,St) 



k=0 



with 



(36) 



and 



(37) 



A, 




Bo 



Bi 




504 
-5832 
2592 




12204 
-13608 





2592 
4428 
-13608 



-13608 
34992 



\ 






-13608 
34992 J 



,B 2 



(0). 



In this affine subspace we want to find a polynomial which satisfies the inequal- 
ities (c) and (d) from Theorem 14. II and which at the same time has a representation 
of the form ([TBI with positive semidefinite matrices F k . Hence, we are left with 
the problem of finding a matrix in the intersection of an affine subspace with the 
cone of positive semidefinite matrices which is a basic task in semidefinite pro- 
gramming. Since this problem is not known to be in NP — in fact it is the major 
open problem in the theory of semidefinite programming — it is not a priori clear 
that a solution exists which one can nicely describe. 

We solved these two semidefinite programming problems separately and we 



used the numerical software csdp [7] for this task: If 0.28 



< 



< 



0.68, then 



F 1 satisfies (c). If 0.18 ^ 7 ^ 0.38, then F 1 has a representation of the form (fTBl 
with positive semidefinite matrices. We make the Ansatz 7=3 and try to find a 
nice representation. For this we solve the semidefinite feasibility problem 



(38) A k + ~B k + PxK ltk + f3 2 K 2 , k + foK 3jk + /3 4 K 4 ,fc >= 0, 



which luckily happens to have the solution (3\ = P3 = P4 



k 

and 2 



0,1,2, 
: 2000. 
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To show uniqueness we first derive the three points distance distribution a of a 
(4, 10, 1/6) spherical code C which is defined by 

(39) a(x, y, z) = T^j l{( c ' c '' c ") £ C 3 : c ■ c = x,c ■ c" = y,c ■ c" = z}\. 

Since —2/3 and 1/6 are the only roots of the polynomial F(x,x, 1) — B, these are 
the only inner products which can occur among distinct points in C. This enables 
us to use (iv) and (v) together with the relations 

a(x, y, z) = a(a(x, y, z)), for all permutations a of x, y, z, 
a(l,l,l) = 1, 

(40) E (w )eD«( a; .!/. z ) = 100 . 
EseI-1,1] a ( x , x A) = 10, 

to determine a by solving a system of linear equations: It is 

a(-2/3, -2/3, 1/6) = 6, a(-2/3, -2/3, 1) = 3, 

(41) a(-2/3, 1/6, 1/6) = 12, a(l/6, 1/6, 1/6) = 18, 
a(l/6, 1/6, 1) = 6, a(l,l,l) = l. 

Now by |[T5l Theorem 5.5] C is a spherical 2-design. By lfl5l Theorem 7.4] it 
carries a 2-class association scheme whose valencies and intersection numbers are 
uniquely determined. In fact it is a strongly regular graph with parameters u = 10, 
k = 3, A = 0, yu = 1. This uniquely defines the Petersen graph which finishes the 
proof of the uniqueness. 
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